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1. INTRODUCTION 



Consider the focusing cubic nonlinear Schrodinger (NLS) equation 

iut + Uxx + '^\u\^u = 0, (1.1) 

where the subscripts denote the appropriate partial derivatives. The NLS equation is 
important for many reasons [1-3,5,31,37]. It arises in many application areas such as 
wave propagation in nonlinear media [37], surface waves on sufficiently deep waters [36], 
and signal propagation in optical fibers [24-26]. It was also the second nonlinear partial 
differential equation (PDE) whose initial value problem was discovered [37] to be solvable 
via the inverse scattering transform (1ST) method. 

In this paper we present a method to construct certain exact solutions to (1.1) that 
are globally analytic on the entire a;t-plane and that decay exponentially as a; — > ±oo 
at each fixed t G R. We derive an explicit formula, namely (4.11), and its equivalents 
(4.12), (5.14), and (6.9), in order to write such solutions in a compact form utilizing 
matrix exponentials. These solutions can alternatively be expressed explicitly as algebraic 
combinations of exponential, trigonometric, and polynomial functions of x and t. We also 
present an explicit formula, namely (5.6), and its equivalents (6.14) and (6.15), for the 
magnitude of such solutions. 

The idea behind our method is similar to that used in [10] to generate exact solutions 
to the Korteweg-de Vries equation on the half line, and we are motivated by the use 
of the 1ST with rational scattering data. This involves representing the corresponding 
scattering data in terms of a matrix realization [11], establishing the separability of the 
kernel of a related Marchenko integral equation by expressing that kernel in terms of a 
matrix exponential, solving the Marchenko integral equation algebraically, and observing 
that the procedure leads to exact solutions to the NLS equation even when the input to 
the Marchenko equation does not necessarily come from any scattering data. 
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For the general use of rational scattering data in inverse scattering theory, the reader 
is referred, for example, to [8,9,17] and the references therein. 

Our method has several advantages: 

(i) It is generalizable to obtain similar explicit formulas for exact solutions to other inte- 
grable nonlinear PDEs where the 1ST involves the use of a Marchenko integral equa- 
tion. For example, a similar method has been used [10] for the half-line Korteweg- 
de Vries equation, and it can be applied to other equations such as the defocusing 
nonlinear Schrodinger equation, the modified Korteweg-de Vries equation, and the 
sine-Gordon equation. 

(ii) It is generalizable to the matrix versions of the aforementioned integrable nonlinear 
PDEs. For example, a similar method has been applied in the second author's Ph.D. 
thesis [20] to the matrix NLS equation in the focusing case with a cubic nonlinearity. 

(iii) As seen from our explicit formula (4.11), our exact solutions are represented in a 
simple and compact form in terms of a square matrix A, a constant row vector C, 
and a constant column vector B, where A appears in a matrix exponential. Such 
matrix exponentials can be "unpacked" in a straightforward way to express our exact 
solutions in terms of exponential, trigonometric, and polynomial functions. Depending 
on the size of A, such unpacked expressions may take many pages to display. Our 
explicit formula and its equivalents allow easy evaluation of such unpacked expressions 
and numerical evaluations on such exact solutions, as evident from the examples in 
available Mathematica notebooks [39] . 

(iv) Our method easily deals with nonsimple bound-state poles and the time evolution of 
the corresponding bound-state norming constants. In the literature, nonsimple bound- 
state poles are usually avoided due to mathematical complications. We refer the 
reader to [32], where nonsimple bound-state poles are investigated and complications 
are encountered. A systematic treatment of nonsimple bound states has recently been 
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given [13]. 

(v) Our method might be generahzable to the case where the matrix A becomes a hnear 
operator on a separable Hilbert space. Such a generahzation on which we are currently 
working would allow us to solve the NLS equation with initial potentials more general 
than those considered in our paper. 

Our method to produce exact solutions to the NLS equation is based on using the 
1ST [1-3,5,31,37]. There are also other methods to obtain solutions to (1.1). Such methods 
include the use of a Darboux transformation [16], the use of a Backlund transformation 
[12,14], the bilinear method of Hirota [28], the use of various other transformations such 
as the Hasimoto transformation [15,27], and various other techniques [6] based on guessing 
the form of a solution and adjusting various parameters. The main idea behind using the 
transformations of Darboux and Backlund is to produce new solutions to (1.1) from pre- 
viously known solutions, and other transformations are used to produce solutions to the 
NLS equation from solutions to other integrable PDEs. The basic idea behind the method 
of Hirota is to represent the solution as a ratio of two functions and to determine those two 
functions by solving some corresponding coupled differential equations. A unified treat- 
ment of Hirota's method, the 1ST, and Backlund transformation to obtain soliton solutions 
with simple and multiple poles for the Sine-Gordon equation was given by Poppe by using 
Fredholm determinants [33]. Other techniques may use an ansatz such as determining 
Q{x,t) and M{x,t) by using u{x,t) = e*®^^'*)M(a;, t) in (1.1). For example, trying 

u{x, t) = e*^'^i^+'=^*+'=3)/(A:4a; + k^t + ke), (1.2) 

where kj are constant real parameters and / is a real-valued smooth function, we get an 
exact solution if we choose k2 = l — ki, k^ = ±1, k^ = ^2ki, and / as the hyperbolic secant. 
One can also use the fact that if U{x, t) is a solution to (1.1), so is e^'^^^~'^^^U{x — 2ct, t) for 
any real constant c. Multiplying a solution by a complex constant of unit amplitude yields 
another solution, and hence such a phase factor can always be omitted from the solution. 
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There are many references in which some exact solutions to (1.1) are presented. For 
example, [38] lists five explicit solutions, one is of the form of (1.2) with a constant /, 
the second and third with / as the hyperbolic secant (those are one-soliton solutions with 
simple poles), the fourth being periodic in x, and the fifth is the n-soliton solution. Another 
solution, which is periodic in x, is [6] 



tt(x, t) — ae 



2ia^t 



2cosh(2a26V2^^t) - ^v'2^^cos(V^a6a;) 



(1.3) 



where a and b are arbitrary real parameters. By letting 6 — > in (1.3) we get the solution 

,2, 3 + 16iaH - IGaH'^ - ^c?x^ 



2ia^t ' 



u(x, t) — ae . „ „ „ 

Another exact solution which is periodic in x is presented [7] in terms of the Jacobi elliptic 
functions. An exact solution to (1.1) is displayed [22] in the form of a specific matrix 
realization and is shown to be valid for t e [0, e) for some small e and x e [0, +oo). In their 
celebrated paper [37] Zakharov and Shabat list the one- and n-soliton solutions as well as a 
one-soliton solution with a double pole, which is obtained from a two-soliton solution with 
simple poles by letting those poles coalesce. In [32] solitons with multiple eigenvalues are 
analyzed and a one-soliton solution with a double pole and a one-soliton solution with a 
triple pole are listed with the help of the symbolic software REDUCE, by stating that "in 
an actual calculation it is very complex to exceed" higher order poles. With our method 
in this paper we show that such solitons with any number of poles and any multiplicities 
can be easily expressed by using an appropriate representation. Let us also add that some 
periodic or almost periodic solutions can be obtained in terms of two hyperelliptic theta 
functions [29,30], and the scattering data for (2.1) can be constructed corresponding to 
certain initial profiles [34,35]. 

In order to appreciate the power of our method, to see why it produces new solutions, 
and to understand why it produces exact solutions that are either impossible or difficult 
to produce by other methods, let us consider the following. When the matrix size is 
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large (imagine A being a 1000 x 1000 matrix) we have an explicit compact formula for 
an exact solution as in (4.11) or its equivalents (4.12), (5.14), and (6.9). By using a 
computer algebra system we can explicitly express such a solution in terms of exponential, 
trigonometric, and polynomial functions of x and t (even though such an expression will 
take thousands of pages to display, we are able to write such an expression thanks to our 
explicit formula). The only explicit formula in the literature comparable to ours is the 
formula for the n-soliton solution without multiplicities. Our own explicit formula yields 
that explicit n-soliton solution without multiplicities in a trivial case; namely, when A is 
a diagonal matrix of distinct entries with positive real parts, as indicated in (7.1). Our 
explicit formula also easily yields the n-soliton solution with arbitrary multiplicities as a 
special case. Dealing with even a single soliton with multiplicities has not been an easy 
task in other methods; for example, the exact solution example presented in [37] for a 
one-soliton solution with a double pole, which is obtained by coalescing two distinct poles 
into one, contains a typographical error, as pointed out in [32]. 

Our paper is organized as follows. In Section 2 we present the preliminaries and outline 
the Marchenko method to solve the inverse scattering problem for the Zakharov-Shabat 
system given in (2.1), summarize the 1ST for the NLS equation, and list in (2.12) the 
time evolution of the norming constants in a compact form [13], which is valid even when 
bound-state poles may have multiplicities greater than one. In Section 3 we consider (2.1) 
with some rational scattering data, which in turn we express in terms of the matrices A, 
B, C given in (3.5)-(3.7), respectively. In Section 4, we derive the explicit formula (4.11) 
for our exact solutions u{x,t) to (1.1) in terms of A, B, C, and we show that such solutions 
have analytic extensions to the entire xt-plane when the real parts of the eigenvalues of 
A are positive. In Section 5 we independently and directly verify that (4.11) is a solution 
to (1.1) as long as the matrix V{x]t) given in (4.7) is invertible, which is assured on the 
entire xt-plane when the real parts of the eigenvalues of A are positive. In Section 5 we 
also show that \u{x^t)\'^ can be expressed in terms of the logarithmic derivative of the 
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determinant of T{x;t). In Section 6 we remove the positivity restriction on the real parts 
of the eigenvalues of A, and we enlarge the class of exact solutions represented by our 
explicit formula (4.11) or its equivalents (4.12), (5.14), and (6.9). Finally, in Section 7 we 
present some examples showing how our explicit formula easily yields exact solutions to 
(1.1) expressed in terms of exponential, trigonometric, and polynomial functions, and we 
also mention the availability of various Mathematica notebooks [39] , in which the user can 
easily modify the input and produce various exact solutions to (1.1) and their animations 
by specifying A, S, C. 

2. PRELIMINARIES 



Consider the Zakharov-Shabat system on the full line 

-iX q{x 





1 









'i 







(2.1) 



—q{x) iX 

where the prime denotes the x-derivative, A is the complex-valued spectral parameter, q is 
a complex- valued integrable potential, and the bar denotes complex conjugation. There are 
two linearly independent vector solutions to (2.1) denoted by il^{X, x) and 0(A, x), which are 
usually known as the Jost solutions and are uniquely obtained by imposing the respective 
asymptotic conditions 





V'(A,x) = 



J— iAa; 





+ o(l), 
+ o(l). 



X 



X 



+00, 



-oo. 



(2.2) 



The transmission coefficient T, the left reflection coefficient L, and the right reflection 
coefficient R are then obtained through the asymptotics 



ipi^x) = 
(f){X,x) = 



,—iXx 



L(A)/T(A) 



7T(A) 



-iXx 



/T(A) 



e*^^i?(A)/T(A) 



+ o(l), 

+ 0(1), 



X 



X 



-oo. 



-l-CXD. 



(2.3) 
(2.4) 
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For further information on these scattering solutions to (2.1) we refer the reader to [1- 
3,5,31,37] and the references therein. 

Besides scattering solutions to (2.1), we have so-called bound-state solutions, which are 
square-integrable solutions to (2.1). They occur at the poles of T in the upper half complex 
plane C+. Let us denote the (distinct) bound-state poles of T by for j = m+1, . . . , m+n, 
and suppose that the multiplicity of the pole at is given by rij. The reason to start 
indexing the bound states with j = m + 1 instead of j = 1 is for notational convenience. It 
is known [1-3,5,31,37] that there is only one linearly independent square-integrable vector 
solution to (2.1) when A = A ^ for j = m + 1, . . . , m + n. Associated with each such Aj , we 
have rij bound-state norming constants Cjs for s = 0, . . . , — 1. 

The inverse scattering problem for (2.1) consists of recovery of q{x) for a; e R from 
an appropriate set of scattering data such as the one consisting of the reflection coefficient 
i?(A) for A e R and the bound-state information {Aj, {cjs}s=o^}j=tn+i- This problem can 
be solved via the Marchenko method as follows [1-3,5,31,37]: 

a) Prom the scattering data {i?(A), {Aj}, {c^s}}, form the Marchenko kernel VL as 

1 /-oo rn+n rij-l 

-^-^ i=m+l s=0 

b) Solve the Marchenko equation 

poo poo 

K(x,y) - Q,{x + y) + / dz ds K{x, s) Q.{s + z) Q,{z + y) = 0, y>x. (2.6) 

Jx J X 

c) Recover the potential q from the solution K{x,y) to the Marchenko equation via 

q{x) = -2K{x,x). (2.7) 

d) Having determined K{x,y), also determine 

pOO 

G{x,y):=- dzK{x,z)n{z + y). (2.8) 

J X 
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Then, obtain the Jost solution 'ip{X,x) to the Zakharov-Shabat system (2.1)-(2.2) via 








K{x,y) 




+ dy 


J X 





(2.9) 



Note that |g(a;)P can be calculated from (2.7) or equivalently by using [37] 



/ dz\q{z)\'^ ^ -2G{x,x), |g(x)p = 2 

J X 



dG{x, x) 
dx 



(2.10) 



The initial-value problem for (1.1) consists of recovery of u{x, t) for t > when u{x^ 0) 
is available. When u{x,Q) = q{x), where q is the potential appearing in (2.1), it is known 
that such an initial- value problem can be solved [1-3,5,31,37] by the method of 1ST as 
indicated in the following diagram: 



, ^ direct scattering r d/ \\ r\ i r n 

U{X,0) > {R{X),{Xj},{Cj,}} 



solution to NLS 

u{x, t) 



time evolution 

{RiX;t),{Xj},{cjsit)}} 



inverse scattering 

The application of the 1ST involves three steps: 

i) Corresponding to the initial potential q{x), obtain the scattering data at t = 0; namely, 
the reflection coefficient -R(A), the bound-state poles Xj of T(A), and the norming 
constants cjs- 

ii) Let the initial scattering data evolve in time. The time-evolved reflection coeflicient 
i?(A; t) is obtained from the reflection coefficient R{X) via 



R{X;t) = R{X) e 



(2.11) 



The bound-state poles Xj and T{X) do not change in time. The time evolution of the 
bound-state norming constants Cjs(t) has been known when s = as 

CjQ{t) = Cjo e^*'^^*, j = n+ 1, . . . jin + n. 
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The time evolution of the remaining terms has recently been analyzed in a systematic 
way [13], and the evolution of Cjs{t) is described by the product of e^*^^* and a 
polynomial in t of order s; we have [13] 

[c,-(,._i)(t) ... c,o(t)] = [c,(,._i) ... Cjo]e-^'^"\ (2.12) 

where Aj is the matrix defined in (3.3). See also [32], where a more complicated 
procedure is given to obtain Cjait). 

iii) Solve the inverse scattering problem for (2.1) with the time-evolved scattering data 
{i?(A; t), {Aj, {cjs(t)}glo^}J!l^_i_i} in order to obtain the time-evolved potential. It 
turns out that the resulting time-evolved potential u{x,t) is a solution to (1.1) and 
reduces to q{x) at t = 0. This inverse problem can be solved by the Marchenko 
method as outlined in Section 4 by replacing the kernel fl{y) with its time-evolved 
version il{y;t), which is obtained by replacing in (2.5) -R(A) by R{X;t) and Cjg by 

Cjs{t). 

3. REPRESENTATION OF THE SCATTERING DATA 

We are interested in obtaining explicit solutions to (1.1) when the reflection coefficient 
i?(A) appearing in (2.4) is a rational function of A with poles occurring in C""". For this 
purpose we will use a method similar to the one developed in [10] and already applied to the 
half-line Korteweg-de Vries equation. We will first represent our scattering data in terms 
of a constant square matrix A, a constant column vector B, and a constant row vector C. 
We will then rewrite the Marchenko kernel fi(j/) given in (2.5) in terms of A, B, C. It will 
turn out that the time-evolved kernel n(y; t) will be related to ^{y) in an easy manner. By 
solving the Marchenko equation (2.6) with the time-evolved kernel ri(j/;/;), we will obtain 
the time-evolved solution K{x^ y; t), from which we will recover the time-evolved potential 
u{x,t) in a manner analogous to (2.7). 

In this section we show how to construct A, B, C from some rational scattering data 
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associated with the Zakharov-Shabat system. We show that our exact solutions can be 
obtained by choosing our triplet A, B, C as in (3.5)-(3.7), where Xj are distinct and 
Cj{nj-i) 7^ for j = 1, . . . , m + n. 

When the rational -R(A) has poles at A-,- in C"*" with multiplicity nj for j = 1, ■ ■ ■ , m, 
since -R(A) ^ as A ^ oo, the partial fraction expansion of R{X) can be written as 



m nj 



(3.1) 



(A- A,)*' 

J = l s = l ^ J' 

for some complex coefficients r^g. Note that we can represent the inner summation in (3.1) 
in the form 



E 



—lyr 



where, for j = 1, . . . , m, we have defined 



J' 



-iXj 










-1 

-iXn — 1 





-1 
zA. 

















-zA. 



- 




"0" 
























-1 







-zA, . 




.1. 



(3.2) 



(3.3) 



so that 



A - iAj = 



{X-iAj)-' = 



rA-A, 











A- A. 























i 

A-A, 

















A- A, 









A-A,J 







(-«) ^ 




(A-A,)2 


(A-A,)3 • 




(A-A,.)"^" 


1 











(A-A,)2 • 
1 


■ (A-A,)"^-^ 


(A-A,.)"^-- 


A-A,- 


■ (A-A,)"^-^ 


(A-A,)"^- 



A-A, 



(A-A,)^ 



A-A, 
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We remark that the row vector Cj contains rij entries, the column vector Bj contains rij 
entries, and Aj is an rij x rij square matrix, {—Aj) is in a Jordan canonical form, and that 
(A — iAj)~^ is an upper triangular Toeplitz matrix. 

As for the bound states, for j = m + 1, . . . , m + n, let us use (3.3) to define the rij x rij 
matrix Aj and the column nj-vector Bj, and let Cj be the row n^- vector defined as 

Cj := [Cj(n,._l) • • • Cjo] , 

so that the summation term in (2.5) is obtained as 

'^J-l s ■ roo 

E ^^>!r = -h / d^Cj{\ - iAj)-'Bje^'y, y > 0. (3.4) 



S=0 



Now let us define the p x p block diagonal matrix A as 





[Ai 


. 


" 







A2 . 





A : = 





. 







. 


. 





where p is the integer given by 



m+n 



rin 



Similarly, let us define the column p-vector B as 



B := 



B2 



B. 



m+n - 



and the row p-vector C as 



(3.5) 



(3.6) 



C:—[Ci C2 ... Cm+n] ■ 



(3.7) 



Without loss of generality we can assume that Aj for j = 1, . . . ,m + n are all distinct; 
in case one of Xj for j = 1, . . . , m coincides with one of Aj for j = m + 1, . . . , m + n, we 
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can simply combine the corresponding blocks in (3.5) to reduce the number of blocks in A 
by one. In case more such Xj coincide, we can proceed in a similar way so that each block 
in (3.5) will be associated with a distinct Xj. Similarly, we can combine the corresponding 
blocks in each of (3.6) and (3.9) so that the sizes of B and C will be compatible with the 
size of A. 

Consider the function P{X) defined as 

P{X) := -iC{X-iA)-^B, XeC, (3.8) 

with the triplet A, B, C, where the constant matrices A, B, C have sizes p x p, p x 1, and 
1 X p, respectively, and the singularities of -P(A) occur at the eigenvalues of iA. Such a 
representation is called minimal [11] if there do not exist constant matrices A, B, C with 
sizes p X p, p X 1, and 1 x p, respectively, such that -P(A) = —iC{X — iA)^^B and p < p. 
There always exists a triplet corresponding to a minimal representation. It is known [11] 
that the realization with the triplet A, B, C is minimal if and only if the two p xp matrices 
defined as 



colp(C, A) 



C 
CA 

CA^-^ 



roWp{A,B):=[B AB ... A^-^B], (3.9) 



both have rank p. 

The following theorem shows that, for the sake of constructing exact solutions to (1.1), 
it is sufficient to consider only the triplet ^, S, C given in (3.5)-(3.7) with distinct Xj for 
j = 1, . . . , m + n because any other triplet A, B, C with sizes p x p, p x 1, and 1 x p, 
respectively, can be equivalently expressed in terms of A, B, C. 

Theorem 3.1 Given any arbitrary triplet A, B, C with sizes p x p, p x 1, and 1 x p, 
respectively, there exists a triplet A, B, C having the form given in (3.5)-(3.7), respectively, 
which yields the same exact solution to (1.1). The construction of A, B, C can be achieved 
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by using 

A = MAM-\ B = MSB, C = CMS, (3.10) 

where M is an invertible matrix whose columns consist of the generalized eigenvectors of 

A, the matrix S is an upper triangular Toeplitz matrix commuting with A, and the complex 
entries of C are chosen as in (3.10). 

PROOF: Since {—A) is in the Jordan canonical form, any given A can be converted to 
A by using A = MAM~^, where M is a matrix whose columns are formed by using the 
generalized eigenvectors of {—A). Next, consider all matrices S commuting with A. Any 
such matrix has the block diagonal form 





-Si 


. 


■ 

































S := 





. 





, Sj : — 








• aj3 


(3.11) 




. 


. 






. 










where nj is the order of the pole Xj for j = 1, . . . , m+n, and the constants ajg are arbitrary. 
We will determine such ajg and hence S itself by using M~^B = SB. Note that SB is 
the column p-vector consisting of m + n column blocks, where the jth block has entries 
aji, . . . , ajrij. Thus, S is unambiguously constructed from M and B. Having constructed 
M and S from A and B, we finally choose the complex entries in the matrix C appearing 
in (3.7) so that C = CMS. Let us now show the equivalence of the representation with the 
triplet A, B, C and that with the triplet A, S, C. From (3.8) we see that we must show 

-iC{\ - iA)-^B = -iC{\ - iA)-^B. (3.12) 

Since SA = AS and MA = AM, we also have 

S{X-iA)-^ = {X-iA)-^S, M{X-iA)-^ = {X-iA)-^M. (3.13) 

Replacing C by CMS on the left hand side of (3.12) and using (3.13), we establish the 
equality in (3.12). Similarly, replacing C by CMS on the right hand side of (4.2) and 
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using MA = AM and SA = AS and (3.13), we prove that ri(j/;t) remains unchanged if 
A, B, C are replaced with A, B, (7, respectively, in (4.2). Hence the triplet A, B, C and 
the triplet A, B, C yield the same solution to (1.1). I 

Note that the invertibility of S is not needed in Theorem 3.1. On the other hand, 
from (3.11) it is seen that S is invertible if and only if ajn^ 7^ for j = 1, . . . , m + n. In the 
rest of this section we will give a characterization for the minimality of the representation 
in (3.8) with the triplet A, B, C given in (3.5)-(3.7). We will show that as long as Xj 
are distinct and Cj(^nj-i) 7^ in (3.7) for j = 1, . . . ,m + n, the triplet A, B, C given in 
(3.5)-(3.7) can be used to recover in the form of (4.11) our exact solutions to (1.1). First, 
we need a result needed in the proof of Theorem 3.3. 

Proposition 3.2 The matrix roWp{A, B) defined in (3.9) is invertible if and only if Xj for 
j = 1, . . . ,m + n appearing in (3.5) are distinct. 

PROOF: It is enough to prove that the rows of roWp(^, B) are linearly independent if and 
only if Aj for j = 1, . . . ,m + n are distinct. We will give the proof by showing that a 
row-echelon equivalent matrix T defined below has linearly independent rows. Using (3.5) 
and (3.6) we get 



roWp(A, B) — 



roWp(^i,Si) 
rowp (^2,-82) 



.rOWp(^^+^, Bm+n) - 

With the help of (3.3) we see that the Uj x p matrix TOWp{Aj, Bj) is given by 

(-1)^"^ 
(-1)^-2 {-ir-\p-l){iX,) 













-1 



2iX^ 



{-ly-^iXjY-^ {-ly-^iXjf-^ 



.1 -iXj {iXjf . 

where we observe the binomial expansion of {—iXj — ly in the (s — l)st column. Put 
a{k) := #{j : nj > k}, i.e. the number of Jordan blocks of A of size at least k. Then, 
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m + n = cr(l) > cr(2) > cr(3) > .... By reordering the rows of roWp(^, B) we obtain a row- 
equivalent pxp echelon matrix T such that T^i = for r > a-(l), Tr2 = for r > (j(l)+cr(2), 
Tr3 = for r > a{l) + a (2) + (j(3), etc., while the submatrices consisting of the elements 

Trs for r = a{l) H h cr(A; — 1) + 1, ... , cr(l) H h cr(A;) and s = k,k+ 1, . . . ,p have the 

form 



(3.14) 



where apart from a sign, the coeflBcients aks are the binomial coeflBcients and hence nonzero, 
and the constants //i, . . . , Hcr{k) correspond to a rearrangement of those of —iXj for which 
Uj > k. Since the matrix given in (3.14) can be written as the product of a Vandermonde 
matrix and a nonsingular diagonal matrix, its rows are linearly independent if and only if 
Xj with Hj > k are distinct. From the echelon structure of the matrix T it then follows 
that all the rows of T, and hence the rows of roWp(^, B) are linearly independent. | 

Theorem 3.3 The triplet A, B, C given in (3.5)-(3.7) corresponds to a minimal repre- 
sentation in (3.8) if and only if Xj are all distinct and Cj(^nj-i) ^ for j = 1, . . . ,m + n. 

PROOF: Note that the matrix S defined in (3.11) commute with A, and we have SA — AS 
and SjAj = AjSj for j = 1, . . . , m + n. Let us use a particular choice for Sj by letting 
aji = Cjo, Oij2 = Cji, . . . , ajrij = Cj(„^_i). Thus, S is invertible if and only if Cj(^rij-i) 
for J = 1, . . . , m + n. Let us define the column p- vector 13 and the row p- vector C via 
B = SB and CS = C. As in the proof of (3.12) in Theorem 3.1 we obtain 

-iC{X - iA)-^B = -iC{X - iA)-^B. 

and hence the representation in (3.8) with the triplet A, S, C is equivalent to that with 
^, S, C. Prom the statement containing (3.9) it then follows that our theorem is proved 
if we can show that roWp(A, B) and colp(C', A) are both invertible if and only if Xj are all 
distinct and Cji^^j-i) 7^ for j = 1, . . . , m + n. Below we will prove that roWp(A, B) and 
colp(C', A) are invertible if and only if roWp(^, B) and S are invertible. Our theorem then 
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follows from Propositfon 3.2 and the fact that S is invertible if and only if Cj(^^._i) ^ for 
j = 1, . . . ,m + n. Since SA — AS and SB — B, from (3.9) we obtain 



SmWp{A, B) — roWp(A, SB) — roWp(A, S), 



and hence roWp(A,i?) is invertible if and only if ToWp[A^B) and 5" are invertible. We 
complete the proof by showing that colp(C', ^) is invertible if and only if roWp(^,S) is 



invertible. Define the rij x rij matrix Jj and the p x p matrix J as 



"0 




1" 

1 



J 



"Ji 
J2 




" 







(3.15) 



1 

1 
















m+n . 



where 1 appears along the trailing diagonal of Jj. Let us use the superscript T to denote 
the matrix transpose. Note that 



It can be verified from (3.5) that JAJ = A^ . Using (3.6), (3.7), and (3.15), since C = CS''^ 



Since J is invertible, our proof is complete. | 

4. EXPLICIT SOLUTIONS TO THE NLS EQUATION 

In the previous section we have constructed A, B, C given in (3.5)-(3.7), respectively, 
from some rational scattering data of the Zakharov-Shabat system. In this section we 
solve the corresponding time-evolved Marchenko equation explicitly for a; > in terms of 
such A, B, C. Such solutions lead to explicit solutions to (1.1) via the formula given in 
(4.11). We then show that such solutions have analytic extensions to the entire a;t-plane if 




we get C = B^J. Thus, we have 




mWp{A^, C^) = rowp(A^, JB) = JroWp(A, B). 
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the real parts of the eigenvalues of A are positive, which is equivalent to having e 0+ 
for j = 1, . . . , m + n in (3.3). We also analyze various properties of the key matrices 
Q{x;t), N{x), and T{x;t) that appear in (4.7)-(4.9) and that are used to construct our 
exact solutions. 

For y > 0, with the help (3.2), (3.4), and a contour integration along the boundary of 
C+, we evaluate the kernel ^{y) defined in (2.5) as 

J](y) = Ce-^^B, y > 0. (4.1) 

Note that (4.1) yields a separable kernel for the Marchenko integral equation in (2.6) 
because from 

n{x + y) = Ce-'^'^e-^yB, 

we see that ^l{x + y) is the Euclidean product of the row p-vector Ce~^^ and the column 
p-vector e~^^5. As a result of this separability we are able to solve the Marchenko integral 
equation (2.6) exactly by algebraic means. 

At this point we discuss the time evolution of the scattering data in more detail. Using 
(2.11) we can express the time-evolved Marchenko integral kernel as 

■, poo m+n rij-l ^ 

^ j=m+l s=0 ^' 

where Cjs{t) satisfies (2.12). This time-evolved kernel is seen to satisfy the first order PDE 

ntiy;t) + 4ii}yy(y;t) = 0, 

provided the integral (1 + A^) |-R(A)| exists. Such PDEs for Marchenko kernels have 

been studied in [4] for a variety of nonlinear evolution equations and in [20] for the matrix 
NLS equation. Here we use (4.1) as an initial condition in solving this PDE and write 

Q{y; t) = Ce-^y-'^'^'^B, y>0. (4.2) 
18 



In other words, ^{y;t) is obtained from by replacing C in (4.1) by Ce"^*"^ *. Let 

us use a dagger to denote the matrix adjoint (complex conjugate and transpose). Since 
ri(j/; t) is a scalar, its complex conjugate is the same as its adjoint and we have 

n{y; t)t = 5te-Aty+4z(At)2t^t. (4.3) 
Comparing with (2.6) we obtain the time-evolved Marchenko integral equation as 

/•OO pOO 

K{x,y;t)-^{x+y;t)'^+ dz ds K{x, s;t) n{s + z;t) n{z+y;ty ^ 0, y > x. (AA) 

J X J X 

Using (4.2) and (4.3) in (4.4), we see that we can look for a solution in the form 

K{x, y; t) = H{x; t) e'^^v+^^^^^^^C^ (4.5) 
where H{x;t) is to be determined. Using (4.5) in (4.4), we obtain 

H{x; t) r{x; t) = S+e"^^^, (4.6) 

where we have defined 

T{x;t) :=I + Q{x;t)N{x), (4.7) 
with / denoting the p x p identity matrix and 

POO 

Q{x;t):= / dse-^'^+^^(^')'*CtCe-^^-^^^'*, (4.8) 

J X 

/•OO 

N{x):= dze-'^'BB^e-^^^ (4.9) 

J X 

Using (4.6) in (4.5) we can write the solution to (4.4) as 

K{x, y; t) = Ste-^'^r(x; t)-ie-^'2^+4^(^')'*Ct, (4.10) 

provided r{x;t) is invertible. We will prove the invertibility of r{x;t) in Theorem 4.2. In 
analogy to (2.7) we get the time-evolved potential as u{x, t) = —2 K{x, x; t), and hence the 
solution to (1.1) is obtained as 

u{x, t) = -2Ste-^'^r(a;; t)-ie-^'^+4»(^')'*Ct. (4.11) 
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It is possible [19] to write (4.11) as the ratio of two determinants as 



u{x, t) — 



detr(a;;t) ' 



(4.12) 



where the (p + 1) x (p + 1) matrix F[x] t) is given by 



F{x;t) :-- 



2Ste-^^^ 



We end this section by hsting some useful properties of the matrices Q{x;t), N{x), 
and r(a;; t). 

Proposition 4.1 The matrices Q{x;t) and N{x) defined in (4-8) and (4-9), respectively, 
satisfy 



Q{x; t) = e-^^^+^^^^^)'*g(0; 0) e-^^-4*^'^ N{x) = e-^^7V(0) e" 



(4.13) 



and the integrals in (4-8) and (4-9) converge for all a;,t e R as long as all the eigenvalues 
of A have positive real parts. 

PROOF: By replacing s and z with s — x and 2; — a; in (4.8) and (4.9), respectively, we 
obtain (4.13). From (4.8) and (4.9), we then get 

Q(0;0)= / ds[Ce-^']^[Ce-'^% N{0) = / [e'^^S] [e-^^S]t. (4.14) 
Jo Jo 

If e > is chosen such that the real parts of the eigenvalues of A exceed e, then in any 



0{e ^^) as 2; — > +00. Hence, 



matrix norm || • || we have = 0{e~^^) and 

the integrals in (4.14) converge, and as a consequence of (4.13) the integrals in (4.8) and 
(4.9) converge for all x, t e R. | 

The next theorem shows that the matrix T{x;t) defined in (4.7) is invertible for all 
t e R as long as the eigenvalues of A have positive real parts. In fact, in that case r{x;t) 
has a positive determinant for all x,t eR. 
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Theorem 4.2 Assume that the eigenvalues of A have positive real parts. Then, for every 
x,t & Ti we have the following: 

(i) The matrices Q{x; t) and N{x) defined in (4-8) and (4-9), respectively, are positive and 
selfadjoint. Consequently, there exist unique positive selfadjoint matrices Q{x;ty/'^ 
and N{x)^/^ such that Q{x;t) = Q{x;ty/^Q{x;t)^/^ and N{x) = N{x)^/^N{xy/^. 

(ii) The matrix V{x]t) defined in (4-'^) is invertible. 

(iii) The determinant ofT{x;t) is positive. 

PROOF: In our proof let us write Q and N for Q{x;t) and N{x), respectively. The 
positivity and selfadjointness of Q and N are a direct consequence of the fact that each 
of the integrands in (4.8) and (4.9) can be written as the product of a matrix and its 
adjoint; hence [23] we have proved (i). From the Sherman- Morrison- Woodbury formula 
[23] it follows that 

[/ + Ql/2(gl/2^)]-l gl/2[j ^ (gl/2^)gl/2] 

and hence (I + QN) is invertible if and only if {I + Q^^'^NQ^/'^) is invertible; on the other 
hand, the latter can be written as [I+{Q^/^N^/^){Q^/^N^/^y] due to the selfadjointness of 
Q^/'^ and N^^'^, and hence it is invertible, establishing (ii). From the two matrix identities 



/ 

Q^/^N I 



I 



-Q^/^N I 



I 






I 










/ 







I + Q^/^NQ^/\ 




/ 


o" 




'I + QN 0" 








/ 




I 





it follows that I+QN and (I+Q^/^NQ^/^) have the same determinant. Thus, we have (iii) 
as a result of the fact that the determinant of [/ + (QV2^i/2^(gi/2^i/2^t] ig positive. | 

Proposition 4.3 Assume that the eigenvalues of A have positive real parts. Then, for all 
x,t the matrices Q{x;t), N{x), T{x;t) defined in (4.7)-(4.9) satisfy: 

Q^ = -A^Q-QA, N^ = -AN -NA^, Qt = Ai[{A^fQ - QA% (4.15) 

21 



T^=I + NQ, r-iQ = Q(rt)-\ (T^)-^N = NT-'^. (4.16) 

PROOF: We obtain (4.15) from (4.13), or (4.8) and (4.9), through differentiation. Using 
the selfadjointness of Q and N proved in Theorem 4.2, from (4.7) we obtain (4.16). | 

Theorem 4.4 For every x,t EH, the matrices Q{x; t) and N{x) defined in (4-8) and (4-9), 
respectively, are simultaneously invertible for all x,t G R if and only if the realization in 
(4-1) of Q{y) with the triplet A, B, C is minimal and the eigenvalues of A have positive 
real parts. 

PROOF: From (4.13) we see that it is enough to prove that Q(0;0) and N(0) defined in 
(4.14) are invertible. The integrals in (4.14) are convergent as a result of the positivity of 
the real parts of the eigenvalues of A. If (5(0; 0)g — for some vector g e C^, then from 
(4.14) we see that Ce~"^^g = for all s > 0. By analytic continuation this implies that 
Ce~"^^g = for all s e C and hence 

CA''g = 0, k = 0,l,.... (4.17) 

Similarly, if N{0)h — for some vector h G C^, using (4.14) we conclude that 

B^{A^)''h = 0, k = 0,l,.... (4.18) 

It is known [11] that the realization in (3.8) or (4.1) for the triplet A,B,C is minimal if 
and only if the two matrices given in (3.9) both have rank p, where we recall that the size 
of A is p X p, that of S is p x 1, and that of C is 1 x p. On the other hand, the ranks of 
the two matrices in (3.9) are both p if and only if (4.17) and (4.18) have only the trivial 
solutions ^ = and h = 0, respectively. | 

For any fixed xq e R, by shifting the dummy integration variable in (4.9) we get 

N{x) = e-^(^-^°)iV(xo)e-^^(^-^°), 

and similarly from (4.8) for any xo,to e R we get 

Q{x;t) = e-^'(^-«^o)+4i(At)2(t-to)g(^^.^^)g-A(x-xo)-4iA^(t-to)_ 
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Thus, we have the following observations. 

Corollary 4.5 Assume that the eigenvalues of A have positive real parts. Then, the matrix 
N{x) defined in (4-9) is invertible for all x & R if and only if it is invertible at any one 
particular value of x. Similarly, Q{x; t) defined in (4-8) is invertible for all x,t if and 
only if it is invertible at any one particular point on the xt-plane. 

Proposition 4.6 If the eigenvalues of A have positive real parts, then the matrix V{x\t) 
defined in (4-'7) satisfies T{x;t) I as x ^ +oo. Additionally, ifQ{0;0) and N{0) given 
in (4-14) 0-1^^ invertible, then T{x;t)~^ — > exponentially as x ^ — oo, where I and are 
the p X p unit and zero matrices, respectively. 

PROOF: As stated in Proposition 4.1, since the integrals in (4.8) and (4.9) converge, 
r(x;t) ^ 7 as a: ^ +oo follows from (4.7)-(4.9). To obtain the limit for T{x]t)~^ as 
a; ^ — oo, let us first define 



Y{x;t) :=e^^r(a;;t) e' 



(4.19) 



Using (4.13) in (4.19) we get 



Y{x; t) = Q(0; t) e-2^^iV(0) / + iV(0)-^e^^^Q(0; t)-^e 



-12Axj 



(4.20) 



Note that, from Theorem 4.2 it follows that N{0)-^ and e'^^''Q{0;t)-^e'^^^'' are positive 
selfadjoint matrices. Using the Sherman-Morrison- Woodbury formula [23] as in the proof 
of Theorem 4.2, we see that the inverse of the matrix in the brackets in (4.20) exists, and 
for all a; e R we have 



Y{x;t)-^= / + iV(0)-V^^Q(0;t)-^e 



12Ah 



iV(0)-^e2^^Q(0;t)-^ 



(4.21) 



Further, since the eigenvalues of A and A'^ have strictly positive real parts, for each fixed 
t G R we conclude, as in the proof of Proposition 4.1, that there exists e > such that 



e^^ = 0{e'^^) and 



= 0(6*^^) as a; — > — oo in any matrix norm || • ||. Hence, from 
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(4.21) we see that Y{x; t) -"^ — > exponentially as a; — > — oo, and writing (4.19) in the form 

we also see that r(a;; t)~^ — > exponentially as a; — > — oo. | 

5. FURTHER PROPERTIES OF OUR EXPLICIT SOLUTIONS 

We have obtained certain explicit solutions to (1.1) in the form of (4.11) by starting 
with some rational scattering data for (2.1) and by constructing the corresponding matrices 
A, B, and C given in (3.5)-(3.7), respectively. In this section we will show that (4.11) is a 
solution to (1.1) no matter how the triplet A, S, C is chosen, as long as the matrix T{x; t) 
defined in (4.7) is invertible. For example, from Theorem 4.2 it follows that T{x; t)~^ exists 
on the entire a;t-plane and thus (4.11) is a solution to (1.1) when the eigenvalues of A have 
positive real parts. 

The purpose of this section is threefold. We will first obtain some useful representa- 
tions for \u{x, t)p corresponding to u{x, t) given in (4.11). Next, we will prove that u{x, t) 
given in (4.11) is a solution to (1.1) as long as V[x]t)~^ exists. Then, we will consider 
further properties of such solutions. 

We can evaluate \u{x^t)\'^ from (4.11) directly. Alternatively, we can recover it by 
using the time-evolved analog of (2.10), namely 



J X 



dz\u(z,t)\'' = -2G(x,x;t), \u(x,t)\'' = 2^^^^^^, (5.1) 

ox 



where, in comparison with (2.8), we see that 

poo 

G{x,y;t):=- dzn{y + z;t)^ K{x, z;t)K (5.2) 

J X 

From (4.3), (4.8), (4.10), and (5.2), we get 

G{x, y; t) = -B^e-^^yr{x; ty^Qix; t)e-^^S. (5.3) 
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Using (5.3) in (5.1), with the help of (4.15), (4.16), and 

(r-^), = -r-^r,r-\ (r-^), = -r-^rtr-\ 



(5.4) 



we obtain 

\u{x, t)\^ = 4Ste-^'^r(x; t)-^[A^Q{x; t) + Q{x; t)A] [r{x; t)^]-^e-^^S. (5.5) 

Next we show that \u{x, t)\^ can be expressed in a simple form in terms of the matrix 
r{x;t) defined in (4.7). As indicated in Theorem 4.2, recall that r{x;t) has a positive 
determinant for all x, t e R when the real parts of the eigenvalues of A are positive. 

Theorem 5.1 The absolute square \u{x,t)\'^ of the solution to the NLS equation can be 
written directly in terms of the determinant of the matrix T{x;t) defined in (4-V so that 



PROOF: In terms of a matrix trace, from (5.1) and (5.3) we get 

\u{x,t)\'^ = -2 

where we have used (4.9) and the fact that in evaluating the trace of a product of two 
matrices the order in the product can be changed. With the help of (4.7), (4.15), (4.16), 
and the trace properties we obtain 

tr [T-^QN^] = tr [-A - A^ + {F^y^A + T'^A^ , (5.8) 

tr [T-^Q^N] = tr [-A - A^ + {T^y^A + T^^A^ . (5.9) 
Thus, from (5.7)-(5.9) with the help of (4.7) we get 

2tr [T-^QNa,] = tr [V-'^Q^N + T-^QN^] = tr [p-^r^] , 
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9det r(x; t)/dx 
det r(x;t) 



dx' 



[log(detr(a;;t))] 



(5.6) 



■Ax 



B 



2tr [r-iQiV,]^ 



(5.7) 



and hence 

which can also be written as (5.6), as indicated in Theorem 7.3 on p. 38 of [18]. | 

We remark that (5.6) is a generalization of the formula given at the end of Section 3 
of [37], where the formula was obtained for the n-soliton solution with simple poles. Thus, 
our method handles the bound states with nonsimple poles easily even though nonsimple 
poles have always caused complications in other methods and have mostly been avoided 
in the literature. 

Let us also remark that (1.1) has infinitely many conserved quantities expressed as 
trace formulas. One such trace formula is given in the following. 

Proposition 5.2 When the eigenvalues of the matrix A have positive real parts, the func- 
tion u{x,t) given in satisfies the trace formula 

dx\u{x,t)\'^ = ii[A + A^] = 2 ^ n^ImfAj], (5.10) 

where Xj and rij are the poles in and the corresponding multiplicities, respectively, as 
in (3.3). 

PROOF: From (5.7) and (5.8) we see that 

/oo 
dx \u{x, t)|2 = tr [-A - ylt + (rt)-ivl + V-^A^] 1^^. 
-oo 

As indicated in Proposition 4.6, we have r(a;;t) — > / as a; — > +oo and r(a;;t)~^ — > as 
X — oo. Thus, we get the first equality in (5.10). Using (3.3) and (3.5), we can write the 
trace of {A-\- A"^) in terms of the multiplicities and imaginary parts of as indicated in 
the second equality in (5.10). | 

Theorem 5.3 The function u{x,t) given in (4-. 11) satisfies (1.1) with any p x p matrix 
A, column p-vector B , and row p-vector C as long as the matrix T{x;t) defined in (4-7) is 
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invertible. In particular, if all eigenvalues of A have positive real parts, then u{x,t) given 
in satisfies (1.1) on the entire xt-plane. 

PROOF: With the help of (4.15), (4.16), and (5.4), through straightforward differentiation 
and after some simphfications, from (4.11) we get 

iut = 8Ste-^'^r-i [{A^f + QA^iVjr-^e-^'^+^^^^'^'^Ct, (5.11) 

Ux = 4Ste-^'^r-i[At - QAiV]r-ie-^'^+^^(^')'*Ct, 
Ua,x = 8Ste-^^^r-^[(At)2 - 2QANr-^QAN + 2A^r-^QAN - 2A^r-^A^ 

(5.12) 

+ 2QANT-'^A^ + Qyl2^^]^-le-^^^+^*^^^)'*C^ 
2uu^u = -16Ste-^^"r-i[(AtQ + QA){T^)-\AN + NA^)]T-^e-^^''+'^'^^^^^^Cl (5.13) 

Using (4.16) and (5.11)-(5.13), and noting that u''' = u, we verify that (1.1) is satisfied. 
Let us note that (5.13) could also be obtained directly by multiplying (4.11) and (5.5). | 

Theorem 5.4 Assume that the eigenvalues of A have positive real parts and that the 
matrices (5(0; 0) and N{0) given in (4-14) osre invertible, or equivalently, assume that the 
representation in (3.8) with the triplet A, B, C is minimal and the eigenvalues of A have 
positive real parts. Then, for each fixed t e R the solution u{x,t) given in (4-11) vanishes 
exponentially as x ^ ±oo. 

PROOF: From (4.11) and the fact that T(x; t) ^ I as x ^ +oo, it follows that u{x, t) ^ 
exponentially as x — > +oo for each fixed t G R. Let us write (4.11) as 

u{x, t) = -2B^Y{x; t)-h''''^^'^''C\ (5.14) 

where Y{x;t) is the matrix defined in (4.19). In the proof of Proposition 4.6, we have 
shown that Y{x;t)~^ —>■ exponentially as a; ^ — oo. Hence, from (5.14) we can conclude 
that for each fixed t e R we have u{x,t) — > exponentially as a; — > — oo. | 

Let us remark that, if the eigenvalues of A have positive real parts, when extended 
to the entire a;-axis the solutions given in (4.11) become multisoliton solutions, where the 
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number of solitons, multiplicity of the corresponding poles, and norming constants can 
be chosen at will. This can also be seen by analytically continuing the time-evolved Jost 
solution 'ip{X,x;t) to the entire a;-axis, by using (2.3), (2.9), and 

hm r dyK{x,y-t)e^^^y--\ (5.15) 

by evaluating the integral with help of (4.10), and by observing that the limit in (5.15) 
vanishes. 

6. GENERALIZATION 

In some parts of Sections 3-5 we have assumed that Xj values appearing in (3.3) and in 
the matrix A given in (3.5) are all located in C"*". In this section we relax that restriction and 
allow some or all Xj to be located in the lower half complex plane C~. Our only restriction 
will be that no will be real and no two distinct Xj will be symmetrically located with 
respect to the real axis in the complex plane. This restriction is mathematically equivalent 
to the disjointness of the sets {Aj}™!^" and {Aj}™!^". Under this restriction we will show 
that u{x, t) given in (4.11) is a solution to (1.1) in any region on the xt-plane in which the 
matrix T{x;t) defined in (4.7) is invertible. The only change we need is that Q{x;t) and 
N{x) will no longer be defined as in (4.8) and (4.9), but instead they will be given as in 
(4.13), where we now let 

(5(0; 0) = — / dX{X + iA^)-'^C^C{X - iA)-\ (6.1) 
27r 

N{0) = — [ dX{X-iA)-^BB\X + iA^)-\ (6.2) 

with 7 being any positively oriented simple closed contour enclosing all A-,- in such a way 
that all Xj lie outside 7. 

As the following proposition shows, the quantities given in (6.1) and (6.2) are the 
unique (selfadjoint) solutions to the respective Lyapunov equations 

Q{0;0)A + A^Q{0;0) = C^C, (6.3) 
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AN{0) + N{0)A'^ = BB^. (6.4) 
We note that, using (4.13), we could also write (6.3) and (6.4) in the equivalent form 

Q{x; t)A + A^Q{x; t) = e-^'-+4^(^')'*ctCe-^"-4^^'*, (6.5) 

AiV(a;) + iV(x) At = g-^^SS+e"^^^. (6.6) 

Proposition 6.1 Assume that none of the eigenvalues of A are purely imaginary and 
that no two eigenvalues of A are symmetrically located with respect to the imaginary axis. 
Equivalently, assume that and {Aj}^^"' are disjoint, where Xj are the complex 

constants appearing in (3.3) and (3.5). We then have the following: 

(i) The matrix equations given in (6.3) and (6.4) are each uniquely solvable. 

(ii) The unique solutions Q{0; 0) and N{0) are selfadjoint matrices. 

(iii) The unique solutions are given by (6.1) and (6.2), respectively. 

PROOF: Note that (i) and (iii) directly follow from Theorem 4.1 in Section 1.4 of [21]. 
It is straightforward to show that the adjoint of any solution to (6.3) or (6.4) is also a 
solution to the same equation, and hence the unique solutions Q(0;0) and A'"(0) must be 
selfadjoint. | 

Next, without requiring that all \j appearing in (3.5) be located in C""", we will prove 
that the matrix u{x,t) given in (4.11) is a solution to (1.1) as long as r(x;t) defined in 
(4.7) is invertible. First, we will write (4.11) in a slightly different but equivalent form. 
Define 

A(a;; t) := / + P(x; t)^ Q{0; 0) P{x; t) 7V(0), P{x; t) := g-^^^-^^^'^ (6.7) 

Note that T{x;t) is invertible if and only if A{x;t) is invertible because, by using (4.7), 
(4.13), and (6.7), we see that 

r(a;;t) = e^^^A(a;;t)e-^^^ (6.8) 
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With the help of (6.8) we can write (4.11) in the equivalent form 

u{x, t) = -2B^A{x; t)-^P{x; t^C^. (6.9) 

Theorem 6.2 Assume that none of the eigenvalues of the matrix A in (3.5) are purely 
imaginary and that no two eigenvalues of A are symmetrically located with respect to the 
imaginary axis. Equivalently, assume that {Aj}J!r^"' and {Aj}™!^" are disjoint, where Xj 
are the complex constants appearing in (3.3) and (3.5). Then, the quantity u{x,t) given 
in (4.. 11), or equivalently in any of (4.12), (5.14), and (6.9) is a solution to (1.1) in any 
region of the xt-plane where the matrix A{x;t) defined in (6.7) or equivalently the matrix 
r{x;t) given in (4-7) is invertible. 

PROOF: In our proof let us write u, A, P, Q, N for u{x, t), A{x; t), P{x; t), Q{0; 0), iV(0), 
respectively. Without explicitly mentioning it, we will use the selfadjointness = Q and 
= N established in Proposition 6.1 as well as the fact that P is invertible. Proceeding 
as in the proof of Theorem 5.3, using straightforward differentiation on (6.9) and after 
some simplification we obtain 

iut = 8B^A-^[{A^f + P^QA^PN]A-^P^C\ (6.10) 

where we have used the fact that 

A = / + P^QPN, {A-^)t = -A-^AtA-\ Pt = -AiA^P, AP = PA. (6.11) 

Similarly, by using (6.11) and 

P, = -2AP, {A-% = -A-'A,A-\ 

after some simplifications we obtain 



= 4B^A-^[A^ - P^QAPN]A-^P^C\ 
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(6.12) 



(6.13) 



u^^ = 8SU-i [(ylt)2 _ 2yltA-iylt + p^QA^PN + 2A^K-^P^QAPN 

-^2P^QAPNK-^A^ - 2P^QAPNA-^P^QAPN] A'^P^C^. 
Next, with the help of (6.3) and (6.4) and using \u\'^u = uu^u^ we obtain 
2|w|2-u= -16StA-i [P^QAP{A^)-^AN + P^QAP{A^)-^NA^ 

+P^A^QP{A^)-^AN + P^A^QP{A^)-^NA^ A'^P^Cl 
We see that (1.1) is satisfied, which is verified by adding (6.10), (6.12), and (6.13) side by 

side and by using 

QPN = {P^)-\A- I), {A^)-^N ^ NA-\ NP^Q = {A^ - I)p-\ 

which directly follows from (6.7) and the selfadjointness of Q and N. | 

As the next theorem shows, if we remove the restriction Xj e C+ then the result in 
Theorem 5.1 still remains valid in any region in the xt-plane where T{x;t) or equivalently 
A{x;t) is invertible. 

Theorem 6.3 Assume that none of the eigenvalues of the matrix A in (3.5) are purely 

imaginary and that no two eigenvalues of A are symmetrically located with respect to the 
imaginary axis. Equivalently, assume that {Aj}^^" and {\j}Y=i^ ^'"^ disjoint, where Xj 
are the complex constants appearing in (3.3) and (3.5). Then, in any region of the xt-plane 
where the matrix A{x;t) defined in (6.7) or equivalently the matrix r{x;t) given in (4-7) 
is invertible, the solution u{x,t) given in (4-11) or equivalently in (6.9) satisfies (5.6) or 
equivalently 

\u{x, t)|^ = tr 



\u{x, t)|^ = tr 



d_ 
dx 

d_ 
dx 



ddeiV{x;t)/dx 
detr(a;;t) 

ddeiA{x]t)/dx 



(6.14) 
(6.15) 



det A{x; t) 

PROOF: Let us write u. A, P, Q, N for u{x, t), A{x; t), P{x; t), Q{x; t), N{x), respectively. 
Using the fact that, in evaluating the trace of a product of two matrices we can change 
the order in the matrix product, from (6.8) we obtain 

tr [r-^r,] = tr [A-^A,] , 
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and hence it is sufficient to prove only (6.14). From (4.13) it follows that (6.5) and (6.6) 
are equivalent to the first two equations, respectively, in (4.15). Note that (4.16) is still 
valid and is a direct consequence of (4.7) and the selfadjointness of Q and N. Proceeding 
as in the proof of Theorem 5.1, with the help of (4.15), (4.16), and (5.4) we obtain 



where we have also used (6.5) and (6.6). Using (4.16) and the aforementioned property of 
the matrix trace, we can simplify the right hand side of (6.17) and show that it is equal to 
the right hand side of (6.16). Finally, as indicated in the proof of Theorem 5.1, the second 
equalities in (6.14) and (6.15) follow from Theorem 7.3 on p. 38 of [18]. | 

7. EXAMPLES 

Specific examples of our exact solutions can be obtained from the explicit formula 
(4.11), or equivalently from any one of (4.12), (5.14), and (6.9), by specifying A, B, and C, 
where T{x-^t) is the matrix defined in (4.7). We have made available various Mathematica 
notebooks [39] in which the user can easily perform the following steps and display the 
corresponding exact solution u{x, t) explicitly in terms of exponential, trigonometric, and 
polynomial functions, verify that the resulting u{x,t) satisfies (1.1), and animate |'u(a;,t)|. 

(i) Input the matrices A, B, C. 

(ii) Evaluate the matrix r{x;t) as in (4.7), where Q{x]t) and N{x) are the matrices 
appearing in (4.13). In case all the eigenvalues of A lie in the right half complex 
plane, evaluate (5(0; 0) and A'"(0) explicitly as in (4.14) with the help of MatrixExp, 



tr [T-^r^] = 2tr [-A - A^ + {T^y^A + T'^A^ , 

tr [r-^r^]^ = 4tr [T-\A^f + {T^y^A'^ - T-^A^T-^A^ 

-{T^)-^A{T^)-^A + 2T-^QANT-^A^ 
On the other hand, using the fact that |wp = uu\ from (4.11) we obtain 



(6.16) 



= 4tr [{AN + NA^)r-\QA + A^Q){r^)-^] , 



(6.17) 
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which is used to evaluate matrix exponentials in Mathematica. In case some or all 
eigenvalues of A lie in the left half complex plane, use (6.1) and (6.2) instead in order 
to evaluate explicitly Q{0; 0) and A'"(0), respectively. 

(iii) Having obtained T{x;t), use (4.11) or one of its equivalents (4.12), (5.14), and (6.9) 
to display u{x, t) explicitly in terms of exponential, trigonometric, and polynomial 
functions. 

(iv) Using (5.6) or (4.11), evaluate \u{x,t)\'^ exactly and animate \u{x,t)\. 

(v) As an option, evaluate the quantities iut, Uxx, and 2|'up'u, and verify directly that 
(1.1) is satisfied. 

Example 7.1 The well-known "n-soliton" to (1.1) is obtained when R{\) = and T(A) 
has n simple bound-state poles in C"*". In this case, from (3.5)-(3.7) we see that A, B, and 

C are given by 

-iXi ... 
-1X2 ... 



A = 











-iXn 





"1" 




1 








.1. 



C = [ Ci C2 ... Cn 



(7.1) 



where are distinct and all lie in C^. Using (4.7)-(4.9), the (a, /?)-entries of the matrices 

Q{x;t), N{x), and r{x;t) are easily evaluated as 



N. 



i{\a—\l3)x 



a/3 



Xn — X 



lCaCj3 e 



;(A^j-A„)a:+4i(A|-A„)t 



a — ^13 

n 

Ta/J = Sal3 - ^ 
7=1 



CaC^ gi(2A^-Ac-A;3)a;+4i(A^-A^)t 



(A^ — AQ,)(A-y — A/3) 

where 601/3 is the Kronecker delta. A Mathematica notebook [39] is available, where the 
user can specify n and {A^, Cj}"^i and display the corresponding u{x, t) explicitly in terms 
of exponential, trigonometric, and polynomial functions and animate |'u(a;,t)|. 



Example 7.2 Choosing 

A: 



2 
-1 



B 



C 
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we evaluate Q(0;0) and A'"(0) using (6.3) and (6.4), respectively. Then, with the help of 
(6.7) and (6.9) we obtain 

8e^**(9e-^^ + 16e^^) - 32e^6^*(4e-2^ + Qe^^) 



u{x, t) 



(7.2) 



-128 cos(12t) + 4e-6a; + 166^^ + Sle-^^^ + 64e2a; ' 

Note that one of the eigenvalues of A in this example is negative and the solution in (7.2) 
is not a soliton solution. A Mathematica notebook containing the animation of (7.2) is 
available [39]. 

Example 7.3 Choosing 



2-i 




-1 
2-i 



B 



C=[l + 2i -l + 4z; 



we get u{x,t) = num(a;, t)/den(a;, t), where 

num(a;,t) :=i024e^^^+^*)-2*(^-^*) [(12 - 9i) + lOOt + (5 - 10i)x] 

+ I3i072e^2(x+4t)-2i(x-6t) ^ ^ ^34 + 32z)t - (2 - Ai)x] , 

den(a;, t) :=25 + 65536e^^(^*+^) 

+ 512e^(^*+^) [I2800t^ + 64(20a; + 43)t + 160a;^ + 304a; + 207] . 
The solution in this example can be described as a soliton of double multiplicity, and its 

Mathematica animation is available [39]. 
Example 7.4 Choosing 

1-1 ] [0" 

C=[l 0], 

we easily obtain u{x,t) = num(a;,t)/den(a;,t), where 

num(a;,t) := 32e"2(^"^**) {[-32768a;2 + 524288^^ + 262144ita; - 65536it] 
+ e-^^ [90112t2 + 15872a;^ + 131072t^a;^ + 4096a;^ + 196608a;t^ 
+ 12288x2 + 9216x + 1344 + 1048576^^^ - 32768^^^ - 35840it 
- 61440zta;] + e"^^ [l28t'^ - 8x^ - 24a; - 15 - 112it - GAitx] } , 

34 





"1 


-1 


" 




"0" 


A = 





1 


-1 















1 




1 



den(a;, t) := 262144 + e"^^ [262144a;^ + 589824a;2 + 393216a; + 524288a;3 

+ 67108864^4 + 8388608x^^2 ^ ^22880] + e"^^ [16384^2 + 4096a;4 
+ 1048576t4 + 15360a; + 344064t2 + 24576a;^ + 131072a;^t2 
+ 393216a;t2 + 3648] + e~^^^. 

The solution in this example can be described as a soliton of triple multiplicity. A Math- 
ematica notebook [39] is available for this example and the corresponding animation. 

In Figures 7.1-7.3 we present some snapshots of 1^(2;, t) \ appearing in Examples 7.2-7.4. 
Further examples of exact solutions to (1.1) expressed in terms of exponential, trigonomet- 
ric, and polynomial functions as well as their animations can be obtained with the help of 
available Mathematica notebooks [39]. It can be directly verified that u{x,t) given in the 
above examples all satisfy (1.1). When the matrix size for the A, B, C becomes large, such 
expressions become lengthy and yet can easily be displayed with the help of Mathematica 
or any other symbolic software. 
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Fig. 7.1 Snapshots of \u{x,t)\ of Example 7.2 at t = 0.0, 0.1, 0.2, 0.3, 0.4, and 0.5. 
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Fig. 7.2 Snapshots of \u{x,t) \ of Example 7.3 at t = -0.5, -0.2, -0.1, 0.0, 0.1, and 0.2. 
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Fig. 7.3 Snapshots of \u{x,t)\ of Example 7.4 at t = 0.0, 0.1, 0.2, 0.3, 0.4, and 0.5. 
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